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Abstract 


The possible response of the atmosphere, as simula- 
ted by the two-level Kintz-Aralcawa global general circulation 
model, to a transient (one month) North Pacific sea surface 
temperature (SST) anomaly Is Investigated in terms of the 
energetics both In the spatial and wave-number domains. In 
the first month, the increased heating, evaporation and con- 
densation over the local anomaly region appear to cause an 
increase In the dlabatic generation of eddy available potential 
energy in northern middle latitudes, and the increased tem- 
perature gradient, especially on the northern side of the 
warm pool, was found to Increase the eddy conversion of the 
intermediate scales In northern middle latitudes. In the 
second month, on the other hand, the removal of the anomalous 
SST pattexui appears to cause a reduction in eddy activity in 
northern middle latitudes, with the reduction taking place 
at ell wave numbers. A large response to the transient North 
Pacific warm pool was also fotmd in the southern middle 
latitudes during the second and third months. In general, the 
results indicate that the transient SST variations of reasonable 
magnitude in the North Pacific Ocean can induce a disturbing 
effect on the global energetics both in the spatial and v;ave- 
number domains. 

The ability of the two-level Nintz-Arahawa model to 
simulate the atmospheric energetics is also examined in this 
study. Ej:cept in tiie tropics, the ti-jO-ievel Mintz--Are.kawa iuCdc?l 


vli 



exhibits a reasonable and realistic energy budget* One najor 
defect of the model simulation is caused by the parameteriza- 
tion of convection, which produces too much convective preci- 
pitation in the tropics. The excessive convective precipitation 
in turn affects all the atmospheric energy transformation pro- 
cesses in the tropics. 
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X . Introd ucblon 

The prediction of weather anomalies and climatic 
fluctuations is a basic goal of meteorology. Several hypotheses 
have been proposed to account for anomalous weather and climate. 
For example, irregular variations of solar activity, variations 
in the proportion of atmospheric carbon dioxide, variations 
in the atmospheric turbidity resulting from intense volcanic 
eruptions, and abnormal surface boundary conditions such as 
sea ice, snow cover and sea surface temperature have all been 
proposed as Important factors in weather and climatic varia- 
tions. some possible effects of sea surface temperature 
anomalies in the Pacific Ocean on the subsequent weather 


pntternR part-.l cnT nrly studied >>7 NamlfiR 


( 1963 * 1969 | 1971 ). 


Bjerlcnes (1966,1969) has also Investigated the possible remote 
effects of anomalous sea surface temperature in the equatorial 


Pacific. 


Pecentlyp several conip'atationally stable general 
circulation roclels have been developed (e,g., hint&, 1955» 


Smagorinsky et al,, 1965: Leith, 196?; Kssahara aiid Washington. 
1967; Kolloway and lianabe, 1971). These models provide the 
tools not only for learning the basic mechanisms of the general 
circulation and climatology of the atmosp’.icre, but also for 
studying the possible responses of the ataosphere to various 
factors which may induce variations in weather and climate. 

Some experimoiits concei’iiing the effects of sea su^^e^.t. tcra- 

Spar (i973,a,b,o) 


peiMture (3ST) anomalies were carried out by 
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with the two-level Mlntz-Arakawa global general circulation 
model (Gates et al., 1971)* Certain hypotheticail ano:::alous 
patterns were superimposed on the Pacific SST field in those 
experiments. The imposed anomaly patterns persist either over 
a season or over one month. The responses of the model atmo- 
sphere to SST anomalies were discussed In terms of the monthly 
mean sea level pressure patterns and certain regional indices. 

The results showed that both transient and persistent SST 
anomalies produced significant re^onses in the hemisphere of 
the anomaly as well as In the opposite hemisphere. 

The energetics of the atmosphere have been investi- 
gated Intensively during the past few years. The purpose of 
studying the atmospheric energy budget is- to reveal the basic 
mechanisms for maintaining the large-scale atmospheric circu- 
lation system. The basic nature of the atmospheric energetics 
may be changed due to climatic variations. For example, 
diagnostic studies of the kinetic energy exchange between 
zonal flow and. large-scale eddies for the troposphere north 
of 20‘^N indlcato the existence of on unusual type of circulation 
in January 19^3 (Viln-Ilielsen, Broim and Drake ^ 1964; Murakami 
and Tomntsu, 1965; Brown, 196?). Studies for that month indi- 
cate that kinetic energy ’was transformed from zc:aal to eddy 
form during this winter period, whereas thie exchange mechanism 
is usually in the opposite direction (e.g,, Coi't, 19^4; Saltsman, 
1970), According to Kamilas (1963) r positive SST nnouiolleB in 
the central end eastern North Pacific Ocean were fairly per- 
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siBtent during the last half of 1962 and early 1963. The 

case of January 1963 seems to suggest that SST anomalies may 

1 

produce significant effects on the atmospheric energetics. 

The fact that the sea level pressure fields exhibited world- 
wide responses to the SST anomalies (Spar,l973aj 1973c) seems 
to indicate that the SST anomalies might induce world-wide 
effects on the atmospheric energetics. 

The purpose of this study is to investigate the 
response of the atmosphere to certain SST anomalies from the 
energetics point of view. The case of the transient (one month) 
SST anomaly (Spar, 1973c) was chosen for detailed analysis. 

In this experiment, the anomalous warm pool was located in 
ths North Pacific Ocean, in an area bounded by 22*^ N - kZ° 
and l^O^W - 180* W, and persisted over a period of only one 
month. In the first month, tvro effects of the SST anomaly 
are expected to affect the atmospheric energetics. One is 
the augmented heating, evaporation, and condensation over 
the local anomaly region, which is expected to increase the, 
diabatlc generation of eddy available potential energy. The 
other is the Increased temperature gradient, especially on 
ths northern side of the v^arm pool. According to the baro- 
clinic instability theory (e.g., Charney, 19^7 ) , the increased 
tejT^peraturc gradient is favorable for baroolinic instability. 


1 

Hov,’evcr, it must be noted that we do not expect the sii::i- 
lar effects of Javmary 1963 to appear in t'-e analysis of this 
study. This is because that ths pabtern and duruation of S3T 
aru::;;al 7 and the initial coniltion of tl;a transient SST anomaly 
experiment (Spar, l973c) are different from these of January 
1963 . 
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Slnoe the conversion of eddy available potential energy to 
eddy kinetic energy is responsible for the growth of the dis- 
turbance in a baroclinic atmosphere (e.g.. Brown, 19^9), this 
conversion Is also expected to increase, especially for the 
most unstable waves. On the other hand, in the second month 
the removal of the SST anomaly Is expected to reduce conden- 
sation and produce a cooling effect over the warm pool region. 
In this study, the propagation of these effects and the in- 
fluence on the other energy transformation processes are in- 
vestigated, The results are discussed in section 7, 

In order to evaluate the signif iconce of such a simu 
lation of climate modification due to SST anomalies, it is 
necessary to assess the ability of the modeT to simulate the 
real atmosphere, A comparison of the energetics of the model 
"control run", which vms designed to simulate the "ncrmal" 
atmosphere, with no SST anomaly, with the "observed" atmoi. 
sphere is presented in section 6. 

The energy equatiens integrated over the entire 


global atmosphere have been derived by Lorens (1955) 
Salt2man(i957) • Since then atmospheric energetics have been 
computed exte^isively based on these equations. Because this 
study is concerned primarily with the possible effects cf SST 
anomalic-s on different parts of the atmosphere, spectral 
energy equations for latitude belts instead of the global 


integrated energy equations were used, Follcving the approach 
of X.orcuc ( 1955 } rnd Saltzmcm(i957) > the energy equatio^iC for 
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a latitude belt both In the spatial and wave-number domains 
were derived, and are presented in sections 2 and3- The model 
energetics were computed front the daily mean values of the 
Fourier representations of variables defined at isobaric 
surface, and the harmonics were calculated up to wave number 
15, based on Information at 72 points around a latitude circle. 
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2. Formulations in the sp:itial doiiiain 

2. 1. Fundamental equations in the spatial domain 

By adopting tlie hydrostatic approximation, the governing fundamental 
equations (i.e. , the equations of motion, thermodynamic energi’, and continiiitj’; in the 
(^» 0» ?3 t) cooi'dinate system may be written as foliows: 


u V 

U + U- + ““ U ^ +11) 

t a cos $ X a 0 p 


u 5, .x = o 

P V a y a cos 6 X 


u 


t' a cos 0 X a 0 


V, + “V+U)V +(f + 


u tan 0 ^ 1 . 

( I + )u+ - + y = o 

V ay a 0 


( 2 . 1 ) 


( 2 . 2 ) 


§ + 
P 


RT 

P 


o 


(2.3) 



1 

a cos 0 


+ (V cos 0)^) = 0 


(2.4) 




u _ 

T 

a cos 0 X 


+ - T + u) T - 
a 0 p 


R 




uu T - 


A 

c 


= o 


p 


in which 

t = time, 

X = longitude, 

0 = latitude, 
p = pressure, 
a radius of the eai'Ui, 
u a cos 6 d-VdL, eastward wind speed. 


(2.5) 



7 


, V = ad -^/dr, nor ilr.vard wind sioccd, 

jj = dp/dt, vertical p-velocity, 

§ = geopotential, 

T - temperature, 

li = diabatic heating I'ate, 

g = acceleration of gra\lty, 

R = gas constant for air, 

c - specific heat at constant pressure, 

P 

f = Coriolis parameter, 

X = eastward component of frictional force, 
y = northward component of frictional force. 

2.2. Energy equations in the spatial domain 

Ihe total potential energy of the atmosphere, wMch is the sum of poten- 
tial and imemal energies, is the only source for the Idnetic energy of the whole atmos- 
phere. However, it is not a good measure of the amount of energy available for conver- 
sion to kinetic energy. For a horizontally stratified and stable atmosphere, although 
total potential energy is plentiful, none at all is available for conversion to kinetic 
energy. Thus, in his discussion of the atmospheric energy cycle, Lorenz (1955) de- 
fined the available potential energj^ of the atmosphere as tlie difference between the total 
potential energy of die whole atmosphere and the total potential cnergj" wliich would exist 
imder any adiabatic redistribution of mass to yield a horizontal stable sti'atification. 
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This concept was also discussed by :\Iargules (1^03) in Ids famous papci’ conceniing the 
energy* of storms. 

The available potential energy v/as fomiiilated by Lorenz (1D55) for iso- 
baric siu'faces., Based on this formulation, the available potential energy has been com- 
puted extensively in observational studies. Modifications of Lorenz's (1955) definition 
of available potential energy have been proposed (e.g. , Van Mieghem, 1956; Dutton and 
Johnson, 1967). In tliis study, Lorenz's (1955) approach is adopted to make the energe- 
tics of the model atmosphere comparable to observations. 

Let K and P be the kinetic and available potential energies per unit mass, 
respectively^ which are defined as 

K = (u^ + v^)/2 

P = Cp Y T"^/2 , 

w'here ( )" denotes the deviation from area average over an isobaric surface, (~), 
and Y (see Saltznian, 1957) is given as 

c p -1 

V" (^- -V • ■ 
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Let [ ] denote tlic zonal average, defined as 

1 

[b] ^ i;; J bd\ , 

"o 

* 

and let ( ) denote the deviation from the zonal average, [ ]. Then the kinetic and 

available potential energies may be decomposed into zonal and eddy components, as 
follows, 

where 

= ([u]^ + [v]^)/2 

* * *2 *2 
= ( [ul u + [V] V ) + (u + V )/2 , 

P^= CpY(XrV2 

K and K are the zonal and eddy kinetic energy per unit mass, respectively; P and P 
2 L z E 

are the zonal and eddy available potential energy per imitmass, respectively. 

The energy equations per unit mass over a latitude belt, derived from the 
fundamental equations (2.1)-C.5), may be written as 
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■"£> - V " ^ °z > 

" <'^z* ^E> ^ <’^E- " "e " °E • 

^z>t = - <^z’ ^E> - <^z> V " ""z - ' 

t^E^t = ^^E> " <Pz> V - <^E> " «E ^ • 

where 


( 2 . 6 ) 
(2. 7) 
( 2 . 8 } 
(2, 9) 


Fay = f^(k^)^f2(k^) . 

^ 2 ^ z '> = - " f "’' «'}0 - "* 1 ^ 

F(y) = F^av-F^ay) , 


[V] [V* w* ] J* 


9 


.w=-{ 


1 

a cos 0 


([V K^] cos 0)^ 


‘"Vp} > 


^2<^v=-^2<^y 


F(P^) = F^(P^).F 2 (P^) , 


F, (P ) = - i — - — “ ([v P ] cos 0) . -^- Y ([iu P ]/v) \ , 

1 z L a cos 0 7/ ffi ‘ -7^' ‘'ril * 


0 


y<y> - { rir? <!t ^" ^ (m" } 


F0P^)^FjlP^)-F3(y) 


9 



u 


r 1 


= - t ‘ ^ - Y (1^ Pe ’ / Y)p } 






Y * ■* 

-- |[U V ] 


*2 

cos 0 . r w . . T V 1 . , * * 

“ 1 “ “ Hu ] 


+ (0)*v*][v]p - [u*^] [V] ^ 


(p ,pp> 
z E 




< Pe >'' e > 


= - c 


a J * 

^{ r ^0 [ vV ][ T ]^ m " t *]([ t ];; - 


— [T]" 


c P 
P 


R ” 

“ J [^][ T ] , 


R * * 

- -[.) T ] . 


W 


w. 


E 


([v](?i" cos 0)^ M[c«i[$3";p} , 
="{ rir 0 ( tv %*] cos 0 )^+ > 


=-[u] [X] - [V] [y] , 


* * * * 

°E X ] - [V y ] , 


G ^= Y [ h ][ T ] , 


Gj ^= Y [ h ’' T *] , 


•I * 4 : O-p 

CRZ = ([03 P 1 -i- Y [T] [i« T ]) 


P 
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M * * 2R 

CRE - ([ojP ] - C V [T| [•:• T ]) — 

h p c^p 


F{b) represents tiie ilux corAcrgence of b; <b ,b ) rcpx’cscnri iJie convcrsicn from 
b, to ((b_,b ) =- (b ,b )); W and W are the energy exchanges with surround- 
ings througli zonal and eddy pressure interactions, respectively; D and D are the 

Z ill 

frictional dissipations of zonal and eddy kinetic energy, respectively; G and G are 

Z ii 

tlie generations of zonal and eddy available potential energy, respectively, by diabatic 
heating. CRZ and CRE arc die extra terms -which arise from the simplified definition 
of available potential energy. These two terms, which were neglected in the global 
integrated energy equations (e,g. , Saltzman, 1957), are disregarded in the following 
discussion. Note that for the area north of 20®^ N, the neglect of CRZ and CRE will not 
make the model computation inconsistent mth the observations. 


Eqs. (2,6) - (2, 9) are used to compute the energetics for latitude belts. These 

equations integrated over the entire mass of the atmosphere are essentially the same as 

tliose derived by other investigators (e.g. , Lorenz, 1955; Saltzman, 1957; Oort, 1964), 

and can be used to estinxate the energj^ cycle for the entire atmosphere. Note that the 

definitions of {P ,K ) and G in this study are diffei*ent from tlioso generally used 

(see Lorenz, 1955; Saltzman, 1957). For example, (P ,K ) and G were defined by 

z z z 

I' t» n H 

Saltzman (1957) as — [ jj] [T] and y [h] [T] , respectively. The flux conver- 

gence terms (F(K ), F(K ), F(P ) andF(P )) and pressure interaction terms (W and 

Z L Z L Z 

\^^ ) are disrcgai'dcd in the global integrated energv equations. However, these terms 

ij 
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cannot be disregarded in the energj' budget for a latitude belt. The other transformation 
teitns in this paper are essentially' the same as those generally defined. 
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3, Fo-n-uilarior.s in the vave-niuv=bcr (lon^nin 

3. 1 Fmidaiiieiiial ccjuations in the viive-nin-nbcr domain 

Let and Z(i^) be the Fouriei' transform pair defined as follov.-s; 


CO 



Z (n) e 


inX. 


(3.1) 


and 


Z(ti) 



p '■ , -inX .. 

J C(>de dX . 

o 


(3.2) 


The Fourier transform pairs to be considered in this study are listed i.n Table 1. 


Table 1. Fourier transform pairs 


c(^) 

u 

V 

UJ 

5 

T 

h 

X 

y 

Z(n) 

U 

V 

Q 

A 

B 

H 

. X 

Y 


By multipljlng Eqs. 


(2. 5} by {2 t0 


"in X 
e 


and integi'ating 


aroimd a latitude circle, the spectral forms of the fundamental equations can be obtained 
(Saltzman, 1957), The spectral equations of horizontal motion may be wa-itteu as 
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lU(n))^ - 


) i U$n) U ^-m) + — U^pii) V^-ni) 


1. a cos 0 


a 


m 


+ U <m)n(r\-m) - ^ Ufrn) V^i-m) f 

n a ■ J 


(3.3) 


in 

a cos 0 


A(n) + £ Vgi) - X{n) 


> 


and 

m = -® 


+ V^&n)n(n-m) + ^ U^i) U^i-m) j- (3.4) 

- i A (n)-fU(^)- Y^) . 

a 9 


where m andn are the wave numbers around a latitude circle. The spectral 
hydrostatic equation may be wTitten as 


A (n) = - — B(q) . 
P P 


(3.5) 


The spectral equation of contimiity may be vritten as 




in ,, , , 1 . tan (5 

I I ^ i 


; 


a 


(3.G) 
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The spectral tJicrmodyr.a-inic 0" :r^y equation may be written ;is 

CO 

T~ /• ITT' T_ 

(B(n)) =- ) i TTTTTT U^-m) + - B (in) V (n-m) 

L ( d COo V' d V 

m = -^ 

+ B (m)O(n-m) - — E^m)Otn-m) i- 

P C p P ' j 

+ “ H^) , 
c 

P 


<3. 7) 


3.2 Energy equations iii the wave-number dom ain 

If Z(n) and Q(n) are the Foux'ier transforms of ^(\) and q(X), respect- 
ively, the Parseval theorem may be WTitten as 


2n 


2tt 

Q{\) q(X) dX 

V 

o 


03 




L 


Q^) Z (-m) , 


IJl = -CO 


(3.8) 


By means of Parseval Uieorem, zonal averages of kinetic and available 
potential energies can be cj^ressed as 


[K] = |([u]^+ [v]^) + Ivwl^) . 

n^l 


and 


IP] - 


1 


c 

P 


Y [T] 


Y c Y . 

P 

n-1 


9 
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If we define 

Kfii) = !u^)J^ + \v(n) f 
P(n) = Y . 

then [K] and [P] can be expressed as 

CO 

[K] = K + y K(tl) , 

Z 1—1 

n=l 

Wi 

[P] = Pjj + y Ptn> , 
n=l 

where K^i) is the eddy kinetic energy of the nth component per unit mass in a latitude 

belt, and P(n) is the eddy available potential energy of the nth component per unit mass 

» . 00 

in a latitude belt. Note tliat K_ = Y K^), and P_ = V P(n), 

E Z_) 

n=^l n=l 

Saltzman (1957) derived the spectral energy equations for the entire 
atmosphere. However, as mentioned before, in order to be able to see tlie effects 
of SST anomalies on different parts of the atmosphere, the calculations were based on 
the spectral energy equations for a latitude belt. 

After multiplying Eqs, (3.3) and {3,4)byU(-n) and V(-n), resj^ectiveiy, 
and the equations for (U(-n)) and (V (-n) ) by U(n) and V(n), respectively, tlie eddy 
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Idnetic energy eciuation for the nth component can be obtained b 3 ' adding the resulting 
four equations (see Saltzman, 1957). Similarly, multipljlng Eq. (3, 7) by B (-n) and 
the equation for {B (-n) ) by B^), respectively, tlie eddy available potential energy 

. L 

equation for the nth component can be obtained by adding the resulting two equations. 
These spectral energy equations for a latitude belt may be wTitten as 


(K(a)h = F(K(n)) + M (n) + L(h) - C(n) + W(n) + D^) , (3.9) 

(P(n)L= F(PW) + Q(n) + S<n) - C(q) + G(n) , (3.10) 

where 


CO 

r-> 


m = -“ 

7^0 


, , TT \ V "N 

CC ^ *vu ~i~ '^vy 0 


0 


( 'i: (m,n) + . > , 

V 2 'uju 2 ‘ujv y J 


F(P{n)}= - 


I 


r 1 

C Y ' ' ' 

p L a cos 0 




m = 

7^0 


+ 


y B(m) , 

2 ’^(UT 


<m,n)^ j 
P 


> 



M 01} = - 


cos 6 


Qfn) - 

LOi) ^ 

1 ^ 2 ^) = 
S(n) = 

S^(ti) = 


^ liJ . '1 

- 7 ~ i [v]^ 


P) [V] 


tan 0 1 

~ y 

a j 




R 


[T]- 


•N 1 

J/ 


Lj^W + L^ftl) 


wu 

I 


*n„ + 7 t,„. Jn.n) 


m = -w 
^0 


a cos 0 'uu. a 'vn 

A 0 


<JJU 


4n,n) - , V«n)(' -1 


a cos 0 ’uV' 

A 


&m,n) 


+ - (m,n) + * (m,n) + 

UJV a "uu yj 


a 


l,„(m,nA!- , 


= F(KOi)) , 
S^Oi) ^2^^ * 


^ a cos 0 ^ &i^>ii) 


m = _co 

4:0 


a vT 


p 



20 


S^^n) = F(P^)) . 

GW = - I X^^fn) , 


= - {r^ } ’ 

D(n) =-x^ W - X^, W . 


GW = Y X^ W , 


in wliich 


<Ji) = [Z(n)Q (-n) + Z(-n) Q^)] , 
iif^^ (m,n) = [ Z(n-m) Q(-n) + Z(-n-m) Q(a)l , 


mth Z^) and Q^) beir^' the Fourier transforms of C^{\) and q{X), respectively. 
F(K(n) ) is tlie flux convergence of eddy kinetic energy of the nth component; F(P^) ) 
is the flux convergence of eddy available potential energy of the nth component; M (n) 
is the transfer from zonal Idnetic energy to eddy kinetic energy of the nth component; 
Q(n) is the transfer from zonal available potential energy to eddy available potential 
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energy of the nth component; L()i) is the kinetic energy transfer from other wave com- 
ponents to Uie ntli component tlirough nonlinear wave interations; S (n) is the available 
potential energy transfer from other wave components to the nth component through 
nonlinear wave interations; C(n) is the eddy conversion from availaJ:)le potential energy 
to kinetic energy for the nth component; W(n) is the energy exchange ivith the surround- 
ings through eddy pressure interation of the nth component; D^i) is the frictional dissi- 
pation of eddy kinetic energy of the nth component; G(n) is the generation of eddy avail- 
able potential energy of the nth component- 

Eqs. (3, 9) and (3, 10) are used to compute the spectral energetics for latitude 
belts. These two equations together with Eqs. (2. 6) and (2, 8) integrated over the 
entire mass of the atmosphere are essentially the same as those derived by Saltzman 
(1957), and can be used to estimate the energy cycle for the whole atmosphere in the 
wave-number domain. 

The relationships between the energy transformation terms in the spatial 
domain and those in the wave-number domain are given in Table 2. 
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4, Brief descriptions of the p;eneral circulation model and 

the transient SST anomaly experiment 

A transient SST anomaly experiment was carried out 

by Spar (1973c) with the two-level Mlntz-Arakawa global 

general circulation model. The details of the model structure 

are described by Gates et al. (I97i)- A very brief description 

of the model and the experiment is given below. 

The model is a two-layer global tropospheric model, 

is bounded by the earth’s surface and the 200 mb level. 

The primitive equations are solved on a spherical grid of 

5° of longitude by 4° of latitude. The large-scale mountain 

system is Incorporated by using cr-coordinate. ( cr = (p - p«p)/ 

(p - P ). where p is the pressure of the tropopause and is 
o T T 

taken to be 200 mb, and is the pressure of earth's surface). 
Land and sea distributions are also Included in the model. 
Temperatures of land are determined from a surface energy 
balance condition. Sea surface temperatures are prescribed as 
the climatological annual mean values. The model also Includes 
a water cycle (with clouds and precipitation), parameterization 
of convection, radiation (including contimiously varying solar 
distance, declination and zenith angle) (Arakawa et al., 1968), 
and a prescribed seasonal variation in the latitude of the 
northern hemisphere snow line. 

In the transient SST anomaly experiment, the 
anomalous warm pool was chosen, to locate in tlie North Pacific 
Ocean, in an area bounded by 22^N-^4-2^^N and l40*^W‘-l80^V/, 
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An anomalous pattern from 2 C to 6 C vras added to the clima- 
tological mean annual sea surface temperature of this area. 

The anomaly on the perimeter grid points. of this area was 
set at 2^C, increasing to 4^C at points one gridpoint in from 
the perimeter, and to a ma3clniuni of 6^C at points two grid- 
points in from the i>erlmeter. The warm oceanic pool was allowed 
to persist for the first month of the anomaly run, after which 
time the climatological mean annual SST pattern was restored. 
The simulated predictions were carried out for 90 forcast days 
of the winter season. 
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5, Methods of calculations 

History tapes of the experiment with the transient 
SST enomaly (Spar, l973c) obtained from the two-layer Mintz- 
Arakawa global general circulation model were used. The general 
circulation model was designed to simulate the troposphere, 
which was assumed to be bounded between 200 mb and the earth's 
surface. In the computation of energetics, since the p-coordinate 
was used, the lower boundary was assumed to be the iOOO mb 
surface Instead of the earth's surface. The troposphere was 
also divided into two layers with the 600 mb surface being the 
interface. 

Variables on an Isobaric surface were interpolated 
from data of history tapes which were computed In 0 * coordinates. 
Horizontal velocities and frictional forces were interpolated 
linearly in T-space, In the interpolations of temperatures, 
geopotentials and temperature lapse rates, it was assumed, 
as in the genered circulation model, that the potential tern- : 
perature was linear in ^ si>ace. For the vertical p-veloclties , 

a parabolic profile in r-space was adopted. 

The diabetic heating includes solar radiation, 
long-wave radiation, condensation and the surface sensible 
heat exchange. The condensation is due to both large-scale and 
cQivective processes. In order to examine the role of various 
heating components in the generation of available potential 
energy, these heating components were separated. In the gen- 
eral circulation model, the large-scale condensation and the 
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surface sensible heat exchange were assumed to occur only in 
the lower layer and were defined on the 0*=3A surface. These 
two heating components were then assumed to be the heating 
components on the 800 mb surface for computing the energetics. 
The rest of the heating components on the 400 mb and 800 mb 
were interpolated linearly in CT-space. 

As in the general circulation model, horizonal 
velocities and frictional forces were defined at latitudes 
from 88®S to with 4^ interval, and the rest of the vari- 

ables were defined from 90^S to 90°N with 4° interval. The 
Fourier components of the variables on the Isobarlc surfaces 
were computed up to wave number 15* The energetics were then 
computed at latitudes from 86 S to 86 N with 4 interval 
using dally mean values of the Fourier components of the 
variables. In calculating the budget of available potential 
energy, the value of / used was the mean value of the stability 
factor over 90 days . 

The summation of L(n) and that of S(n) over 15 vjave 
numbers are in general not equal to zero. This is because 
of the truncation in wave number, of the finite differences 
employed in calculating the values of derivatives of vari- 
ables, and of the nonzero of the CRE. The nonzero sums were 
assumed to be uniformly distributed over all v^eves. The cor» 
rected L(n) was set to be equal to L(n) defined in section 3 
minus a correction which is the summaticri of L(n) over all 
waves devided by the number of waves. The same approximation 



was also adopted earlier by Saltzoan and Plelsher (i960) and 
Yang (1967). 

The energetics Integrated through the entire depth 
of the model atmosphere were computed. Since the vertical 
p-velocities were assiimed to be zero at 200 mb surface, inte- 
grating the divergences of vertical transport terms ( fcOK_) 

z p 

(^^Kg)p , etc.) yields the vertical transports at the lower 
boundary. In order to Incorporate the topographic effect, 
this lower bo^mdary was set to be the earth's surface rather 
than the 1000 mb surface. Six pressure surfaces, l.e,, 700, 
800, 850 , 900,950 and 1000 mb, were chosen to represent the. 
lower boundary for different grid points in calculating the 
vertical transports of available potential energy and geo- 
potential deviations near the earth's surface. If the 
earth* s elevation at a grid point is higher than the height 
of 775 nib in the NACA standard atmosphere, the lower boundary 
for that grid point was set at 700 mb. If it is lower than 
the standard height of 975 21 b, 1000 rob was chosen to be the 
lower boundary for that grid point. 
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6. Enerfcetics of the nodel atmosphere 

6.1. Energ:etlc5 of the control run In the spatial domain 
for the region north of 
Before Investigating the possible effects of a 
transient SST anomaly on the atmospheric energetics, it is 
desirable to examine how well the model atmosphere simulates 
the real atmospheric processes and energy contents. Since 
the general circulation of the control run was designed to 
simulate normal climate, the energetics of the control run 
are compared with the observations. 

The energetics of the atmosphere have been computed 
Intensively based on observed data. However, most of the 
observational studies are confined to the. Northern Hemisphere, 
especially to the region north of 20°K (see Oort, 1964; 

Krueger et al,, 1965; Saltzman, 1970; etc). Since the ener- 
getics might be quite different over different parts of the 
atmosphere, only the results for the control run north of 
20®N are compared with the observations. 

Fig. 1 shows the ninety-day mean energy budget in 
the spatial domain for the winter control run averaged over 
the area north of 20^N, together with observational results 
from Krueger et al, (1965) and Saltzman (1970). The energy 
estimates by Krueger et al. (1965) shown in Fig. 1 are for 
an average of 5 winter (December- February) seasons, and are 
based on data from 85O mb to 500 mb in the region north of 
20 N, extrapolated to include the entire atmosphere. Eased 
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on a variety of sources, Saltzraan (1970) compiled a spectral 
energy budget for the Northern Hemisphere. His results for 
an average "cold season", October - March, transformed into 
the spatial domain, are also shown in Fig, 1 for comparison. 

In general, the so-called observational energetics 
of the atmosphere have been estimated (e.g,, Oort, 1964; 
Saltzman, 1970) by applying the global Integrated energy 
equations to the data in certain latitude bands within the 
Northern Hemisphere (approximately north of 20°N or 15^N), 

The model energetics, on the other hand, were computed from 
the zonal averaged energy equations for all latitude bands^ 
covering the entire global atmosphere below 200 mb. The 
region north of 20^N is an open domain in the model calcula- 
tions, therefore, the model energetics include the energy 
flux convergence and pressure interaction terms, which were 
neglected In the observational studies. According to Fig.l, 

the convergence of zonal available potential energy, P(P ), 

z 

and the zonal pressure Interaction, W , play en important 

z 

role In the energy budget of this region. 

The zonal available potential energy, P_ , of the 
winter control run appears to be too large compared with the 
observations. One reason for this discrepancy is the differ- 
ence of domains between the model end observational computa- 
tions, Squares of the temperature deviations from the global 
mean temperature were used to compute P^ in the model energe- 
tics, while deviations from the limited area mean temperature 
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(e.g. , north of 20^N) were used in the observational studies. 
Another reason is that the model produces excessively high 
temperatiire In the upper troposphere in the tropics. In the 
January climate simulation, Gates (1972) also found that the 
temperatures in the upper level of the tropics were much 
higher than those observed. This Is because the model pro- 
duces too much convective precipitation in the tropics. 

The zonal kinetic energy, t of the model control 
run Is In good agreement with the observations, and so is the 
eddy available i>otential energy, Pg , On the other hand, the 
eddy kinetic energy, Kg , calculated from the model is too 
small, probably as a result of the low horizontal resolution 
of the model. Both Kanabe et »1, (1970) and Wellck et al, 
(1971) found that the eddy kinetic energy was markedly in- 
creased by increasing the horizontal resolution of the model* 
The zonal conversion indicates a conver- 
sion from zonal kinetic to zonal potential energy in the 
model north of 20°N, in good agreement with the result of 
Krueger et al. (1965). The negative of this conversion 

o 

reflects the dominance of the indirect cell north of 20 N. 

On the other hand, the eddy conversion, {,Pg,Kg)- , from eddy 
available potential to eddy kinetic energy, appears to be 
too large in the model compared with the observations. ROwever, 
this conversion was computed from the vertical motion field, 
which was determined dry adiabatically in the observational 
diagnostics. The observational {PpiKr) s^re probably under- 
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estimated due to the underestimated vertical velocities. 

Since the moist convection was included in the model, the 
value of more realistic in the model com- 

putation than in the =*observational" diagnostics* 

The term indicates a transformation from 

eddy to zonal kinetic energy in the model, in good agreement 
with the observations. However, the transformation (Pz.I’e) . 
from zonal to eddy available potential energy, seems too low 
in the model, probably as a result of the low horizontal re- 
solution of the model* 

The generation of zonal available potential energy, 

(>2 t shows good agreement between the model and observations* 
However, for the generation of eddy available potential energy, 
Gg , there is a disagreement in sign. The eddy available po- 
tential energy is generated by diabatic heating in the model, 
but is destroyed in the observational diagnostics. Both G2 
and Gg compiled by Saltzman (1970) are from studies of 
Wiin-Nielsen and Brovm (196O) and Brown (1964), in which the 
two- parameter baroclinic model was used. The two-parameter 
baroclinic model might underestimate the vertical velocities, 
which in turn caused an underestimation of the creation of 
Pg by condensation heating. This is probably the reason why 
the values of in Saltzman’s study appear to be too large 

and are negative in sign. Furthermore, the values of G and 

z 

Gg were determined as residuals by Krueger et al.(l965). 

Thus, these observational values cannot be considered as 
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reliable standards against which to Judge the niodel. 

The kinetic energy dlesipatlons , and Dg , in 
the model are not in good agreement with those of Saltzman 
(1970). The total dissipation, (0^+ Dg), is 6.43 watt 

in the model compared with Saltzman'a residual estimate of 

•2 

3*03 watt m . However, based on five-.year‘s data, Kung 
(1967) suggested a value of 4.94 vjatt for the winter 
(October - March) total dissipation, which is in better 
agreement with the model. But, as the zonal winds in the 
middle latitudes are much higher in the model than in the 
observations (Gates, 1972), it is likely that the zonal dis- 
sipation, , is somewhat too high in the model. 

6,2. Energetics of the control i*un in the wave-number 
domain for the region north of 20^?; 

Fig. 2 shows the ninety-day mean energy budget in 
the wave-number domain for the winter control run averaged 
over the area north of 20 H. Except for eddy generation terms, 
G(n) , the model energy flow diagram in the wave-number doamin 
is qualitatively In good agreement with the observational 
results of Saltzraan (1970), Note that Fig. 2 includes energy 
flux convergence and pressure interaction terms, which were 
neglected in the observational results of Saltzman (19?0). 

The general picture is described below. 

The available potential energy generated by diabatic 
heating takes place mainly in the zonal form (G^), and flows 
to all the eddies primarily through eddy sensible heat trans- 
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fer measured by Q(n), with the dominant transfer in the larger 
eddies (small wave numbers). The available potential energy 
is further transferred from longer waves to intermediate and 
short waves through non-linear wave interactions, S(n). All 
waves convert their potential energy to kinetic energy through 
barocllnlc processes represented by C(n), and supply kinetic 
energy to the zonal flow through barotropic processes measured 
by K(n). Intermediate waves lose kinetic energy, and most 
longer waves and short waves gain kinetic energy through non- 
linear wave interactions, L(n). 

However, the eddy energies, P(n) and K(n), are' 

generally smaller in the model than in the observations, as are 

the conversions from to P(n), Q(n), Also, in the non-linear 

z 

wave interactions, intermediate waves lose less kinetic energy 
In the model than in the observational diagnostics. Some of 
these shortcomings of the model are undoubtedly due to the 
coarse horizontal resolution, as shown by Kanabe et al. (1970). 

Maximum eddy conversion, C(n), in the model was 
found in the ultra-long waves (near n=2), whereas the barocllnlc 
instability theory and observational studies (e,g,, Saltzman and 
Fleisher, I960, 1961 ) Indicate a maximum in the Intermediate 
scale, around n=6. Most of the barocllnlc instability theories 
were based on dry adiabatic models, and the "observational" 
vertical velocities were computed dry adinbatlcally , in the 
study of the effect of sensible heat exchange on the barocllnic 
instability, Haltiner (196?) found that the heating decreased 
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the instability at Intermediate w?wrelengths .and increased the 
Instability at longer wavelengths, Murakami and Tomatsu (l965) 
further showed that the most unstable wave could occur at 
ultra-long waves if sensible heat exchange was taken Into con- 
sideration, The maximum conversion in the model of Manabe et 
al, (1970) also takes place in the ultra-long wave, as in the 
model. The fact that the eddy conversions of the model at all 
wave numbers are larger than those compiled by Saltzman (1970) 
is probably also due to the use of dry adiabatic vertical 
notion in the observational diagnostics. 

In summery the model control run exhibits a rea- 
sonable and realistic energy budget both in the spatial and 
wave-number domains for the region north of 20°K, 

6.3 Latitudinal distributions of the mcdel energetics 

There are more similarities than differences be- 
tween the energetics of the model control run and those of 
the model anomaly run. The essential features of the model 
energetics are reflected in the similarities. Figs. 3-6 show 
the thirty-day mean latitudinal distributions of energy contents 
for the control and anomaly runs for the second £md third 
months. Although not shovm in figures, the chf-racteristics 
of the model energetics in the first month is similar to that 
in the second and third months. There are tvro maxima of zonal 

kinetic energy, K , in the middle latitudes of both hemi- 

z 

spheres reflecting the jet streamis in the respective latitudes 
(Figs. The eddy kinetic encr^:;/, also has maximum 


























values In the northern and southern middle latitudes corres- 
ponding to the activities of the extratroplcal disturbances 
(Pigs, A minor peak of K appears near the equator 

which shows the activities of tropical disturbances. 

The zonal available potential energy, P , has two 

z 

minima in the middle latitudes of both hemispheres, which are 
probably too close to the equator due to the excessively high 
tropical temperatures (Figs. 4 a, 6A). These two minima could 
be shifted further pole-wards by improving the convective 
adjustment to lower the tropical temperatures. In the north- 
ern hemisphere winter, the coldest temperature occurs at the 
north pole, hence, the highest P, is found there. The eddy 
available potential energy, P- , has maxlrrnm values \n the 
middle latitudes of both hemispheres, reflecting the large 
longitudinal temperature variations In the middle latitudes 

(Figs, 4B,6B), However, the peak of P near the equator is 

E 

probably exaggerated and Is overgenerated by the convective 
condensation. 

Except in the model tropics, the patterns of ener- 
gies In the model atmospheres are qualitatively in agreement 
with the observations (e.g., Vincent, 1969). 

Figs, 7-16 shov/ the thirty-day mean latitudinal 
distributions of energy transformations for the second and 
third months. One noticeable feature is that there exists 
a spike within the tropical region for most of the energy 
transformations . There is good reason to believe that these 
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reflect a defect In the nodel, namely the parameterization of 
convection. The excessive convective condensation in the 
tropics could overgenerate the available potential energy, 
flnd could increase the vertical velocities which in turn In- 
tensify the conversion from potential to kinetic energy. 

These results might further affect other atmospheric pro- 
cesses in the tropics* 

4 

In low latitudes the transformation 
from eddy to zonal available potential energy, and agrees 
qualitatively with that' proposed by Starr and Wallace (196^)* 
However, its magnitude is probably too large in the model 
tropics (Pigs. 9A, i4a)* The eddy dissipation, Dg , in the 
model tropics is also probably too large for balancing with 
the exggerated eddy conversion < Pg,Kg ) (Figs* 8B, 113, 13B, 
i6b)* The eddy pressure interaction, Wg , in the model tro- 
pics even reverse in sign with respect to the results in the 
general circulation model of Manabe et al, (1970a)* 

It is generally agreed that zonal available po- 
tential energy is generated by low-latitude heating and high- 
latitude cooling (e.g*, Vincent, 1969). However, the model 
shows that in low latitudes outside the tropical spike, P 
is destroyed by diabatic cooling (Figs. ?A, 12A). Further 
investigation shows that in low latitudes the major heating 
component destroying P is long-wave radiation, which might 
be too large resulting from the excessively high tropical 
temperatures* In the middle: and high latitudes, the patterns 
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of G In the model are qualitatively in agreement with the 
z 

observational results of Vincent (1969) (Figs, ?A, 12A). In 
high latitudes P is primarily generated by the cooling of 
long-wave radiation. On the other hand, the patterns of G„ 
(Figs, 7B, 12B), in the model are determined mainly by con- 
densation processes. Except in the model tropics, the values 
of G_ are generally small compared with G , 

£• Z 

The patterns of ( ^ exhibit three-cell nerl- 

dlonal circulations in both hemispheres (Figs. 8 a, 13A). 
However, the strength of the Hadley circulations is too 
strong resulting from the overproduction of convective pre- 
cipitation in the tropics. 


The eddy conversion. 


/ r> V \ 


/ f 


has maximum 


*r e 1 «« AO 

m MW W 


in the middle latitudes of both hemispheres reflecting the 
extreme baroclinlc instability in these regions (Pigs. 8B, 
l3B)i The difference of lend-sea distributions in both hemi- 
spheres and the winter conditions in the Northern Hemisphere 
enhance the eddy activities in the northern middle latitudes. 
Thus, the peak values of in the northern middle 

latitudes are larger than those in the southern middle lati- 
tudes , 


The available potential energy is found to transfer 
from zonal to eddy in the middle and high latitudes, with 
the maxima in the middle latitudes (Figs, 9A, i4a). This is 
qualitatively in agreement with the observational studies 
(e.go, Oort, 1964; Vincent, 1969; Saltzmon, 1970), 
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The eddy pressure interaction, , acts as a source 
of eddy kinetic energy In the subtropics £ind as a sink in the 
middle and high latitudes (Figs, lOB, 15B), The results in 
the northern subtropical and middle latitudes are qualita- 
tively in agreement with the results of Manabe et al, {l970a). 
In the general circulation model of Kanabe et al. { 1970a), 
also supplies energy to high latitudes, which is contrary to 
the results In this study. 

The zonal dissipation, D , has ra«oclmuin values in 

£t 

the middle latitudes reflecting the fact that maximiua dissipa- 
tion occurs in the vicinity of the Jet stream (Pigs. IIA, i 6A). 
The eddy dissipation, Dg , also has maximum values in the 
middle latitudes, as expected (Figs* IIB, l6B). 

Except in the tropics, the hlntz-Arakawa two layer 
global general circulation model exhibits a reasonable eind 
realistic energy budget. In the tropics the convective pre- 
cipitation was found to have a dominating Influence on the 
model energetics. Thus, improvements on the parameteriza- 
tion of convection may bring the simulated tropical energe- 
tics closer to the reality. 



55 


7. Effects of a transient SS7 anomaly on the enerr-etlcs 

7,1, Effects on the latitudinal distributions of energetics 
The analyses of the sea level pressure pattern and 
600 mb geopotential field by Spar (1973, a, b, c) showed that 
In the first month the responses to the SST anomaly are mainly 
in the Northern Hemisphere. The ^alysis of the atmospheric 
energetics in this study also revealed that in the first month 
the effects are primarily in the northern middle latitudes. 
Since the responses are smaller in the first month than in 
the second and third months, they will be discussed later 
only in terms of regional analysis. Effects on the latitui^ 
dinal distributions of energetics for the second and third 
months are discussed below. 

(A) Effects in the second month 

In the SST anomaly experiment, a positive SST 
anomaly was introduced to the anomaly run on the first day 
of the first month, and was removed from the anomaly run on 
the first day of the second month. In the first month, with 
the insertion of the SST anomaly, the portion of the atmo- 
sphere above the warm pool tends to gain more sensible heat 
and water vapor from the underlying sea surface; the surface 
air above the warm pool is heated gradually by the warm sea. 

In the second month, with the removal of the SST anomaly, 
this portion of the atmosphere tends to reduce surface sen- 
sible heat flujc and evaporation. Therefore, the transient 
SST anomaly has a heating effect in the first month but has 
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a cooling effect in the second month In this portion of the 
atmosphere. These effects then propagate Into different 
parts of the atmosphere through the atmospheric circulations. 
Since the mean motion of the atmosphere are mainly parallel 
to latitude circles and the cyclone paths are toward north- 
east, it is likely that in the second month the cooling 
effect propagates mostly to the other portions of the north- 
ern middle latitudes rather than to other latitude bands. 

During the second month, as expected, the anomaly 
was found to decrease the amount of condensation and eddy 
generation in the northern middle latitudes (Pig. 7B), It 
also reduces the energy conversions, and ( Pg,Kg ) , 

in the northern middle latitudes (Figs 8A, SB), The fojrmer 
indicates a reduction In the strength of the Ferrel cell, in 
which the vertical velocities In the upward branch of the 
cell are decreased. The latter indicates a decrease in the 
baroclinic instability, which further causes a decrease in 
(Fig, 3B). The less active eddies dissipate less energy in 
the northern middle latitudes (Fig. IIB). The conversion 

P,,P ) in the northern middle latitudes also decreases in 
\ z' E ' 

the anomaly run compared with the control ruji (Fig. 9A), ^ 

This together with the reduction in causes a decrease 

P„ (Fig* 4B). Compared with the control run, the zonal/pres- 
h / 

sure interaction, , in the anomaly run appears to^supply 
less energy to the northern middle latitudes (Pig/ lOA), 

On the other hand, in the anomaly run removed less energy 

hi 
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from the northern middle latitudes .and supplies less energy 

to the northern subtropics, compared with the control run 
(Pig. lOB). 

In summary , the transient SST anomaly was found to 
cause a reduction in the atmospheric processes in the north- 
ern middle latitudes during the second month. Apparently, 
this is the response to the second month's cooling rather 
than to the first month's warming. 

In the second month, the anomaly also affects other 
area outside the northern middle latitudes, especially the 
southern middle latitudes. For example, the anomaly appears 
to cause an Increase in , and < K^.K^) . and a 

decrease In ( southern middle latitudes (Figs. 

3A, 3B, 4b, 9A, 9B). However, the response is less pro- 
nounced and less systematic in the area outside the northern 
middle latitudes, 

(P) Effects in the third month 

As in the second month, the sea surface temperatures 
in the third month are the same for both the anomaly and con- 
trol runs. However, the differences of energetics tstween 
the anomaly and control runs are still large. Haximum dif- 
ferences were found in the tropics. In the tropics, all the 
zonal and eddy energies are decreased in the anomaly run, 
compared with the control run (Figs, 5 , 6 ). Temperatures in 
the upper troposphere were found to be much lower in the 
anomaly rim, with a maxircu-n reduction in the tropics. 


This 
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results in a lower global niean temperature, and hence P is 

z 

markedly decreased in the tropics (also in high latitudes) 

(Pig, 6 a), The anomaly also appears to reduce wind velocities 

in the tropics. This can be seen from Pigs. 5A and 5B that 

both K_ and K in the tropics are reduced in the anomaly run. 
z E 

The energy transformations in the tropics are 
severely distorted in the anomaly rmi compared with the 
control run (Figs. 12 - l6). Compared with the control run, 
the generations, conversions and dissipations in the tropics 
appear to be less active in the anomaly run. The precipita- 
tion and latent heat release in the tropics are due primarily , 

to convective processes. Except at 2^S, the precipitation 

o o 

in the third month in the equatorial region (20 N • 20 S) 

appears to be larger in the anomaly run than In the control 

run. However, at 2®S, where all the energy transformations 

show dramatic reductions, the decrease of precipitation in 

the anomaly run was extremely large (about 80,^). Since the 

convective condensation is the major component for generating 

P and P_ in the tropics, the dramatically altered convective 
z b 

condensation markedly changes the patterns of 0 and in 

Z E 

the tropics in the anomaly r\jn (Figs. 12A, 12B), These 
changes further affect all the other energy transf ormation 
processes In the tropics. Because of the imsatisf actory 
parameterization of convection In the model ( Gates, 1972), 
it is very likely that this tropical effect is a computational 
artifact caused by a model defect rather than a true reflec- 
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tion of nature. 

In the third month the influence of the SST anomsily 
in the northern middle latitudes is hot so pronounced euid 
systematic as in the second month. Nevertheless, some res- 
ponses can still te found. Both K^. and P_ in the northern 

£• E 

middle latitudes appear to he larger in the anomaly run than 
in the control run (Figs, 5B, 6B), As can be seen from the 
pattern of { PziKg^ > the indirect Ferrel cell in the north- 
ern middle latitudes is intensified in the anomaly run (Pig. 
13A). Compared with the control run, it was found that Wg 
in the anomaly run transports more energy out of the northern 
middle latitudes (Fig. 15B), However, the conversion from 
Pg to Pg In the northern middle latitudes is reduced In the 
anomaly run (Pig. i4a). 

For the southern middle latitudes, on the contrary, 
the response to the SST anomaly is more apparent and more 
systematic in the third month than in the second month. The 
transient Northern Pacific warm pool appears to cause a de- 
crease in the eddy activity in the southern middle latitudes 
during the third month. This is shown In the reductions of 
Kg , Pg , ^ Pg,Kg) , ^ southern middle 

latitudes in the anomaly run (Figs, 5B, 6e, 13B, i4a, 16B), 
The pattern of V-. also reveals the reduction of eddy activity 
in the southern middle latitudes. The eddy pressure inter- 
action in the anomaly run removes less energy out of the 
southern middle latitudes, and supplies less energy to the 
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subtropics of the Southern Henjisphere (Fig. 15B), The con- 
version < KgjK^) in lovi latitudes of about 40°S - 40^N 
appears to be increased by the transient SST anomaly. 

Although the anomalous pattern of the sea surface 
temperature field persists for a period of only one month, 
the results for the second and third months suggest that 
the transient SST anomaly can Indeed have a disturbing effect 
on the global energetics. 

7.2. Effects on the regional energetics 

It is well known (Mintz et al., 1972) that the two- 
level Klntz-Arakawa model produces a poor simulation of the 
climate in the polar and tropical regions. It is also a 
fact that (except for the dubious tropical response in the 
third month) the transient North Pacific SST anomaly in the 
experiment produces a maximum response in the middle latitudes. 
Because of these two reasons, two zonal sections representing 
middle latitudes were particularly chosen for detailed analy- 
sis, These two zonal sections are 28^N - 60^N and 28 S - 60 S. 
The effects in these two areas are discussed in tei*ms of 
gross regional energetics both in the spatial and wave-number 
domains . 

(A) Effects in the northern r*iddle latitudes 

The thirty-day mean energetics for each month of 
the control and anomaly runs in the zonal section 28^N - 60°N 
are shown in Table 3- Figs. 17, 18, and 19 illustrate the 
thirty-day meaii spectral distributions of K(n), P(n), K(n), 
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Table 3 « Thirty-day mean energetics for each month, of the 

control and anomaly runs in the zonal section 

28®N - 60^N, The unit of energy content is 10^ 

2 

joules m and that of energy transformation is 
10*^ watt 


Month No. 

Control 

1 

Anomaly 

1 

Control 

2 

Anomaly 

2 

Control 

3 

Anomaly 

3 

Pz 

418 

422 

468 

517 

354 

339 


122 

131 

138 

129 

89 

104 

K 

z 

169 

171 

183 

192 

147 

. i49 

’'e 

100 

100 

113 

106 

73 

82 

Gz 

177 

152 

295 

321 

230 

177 


77 

108 

89 

0 

84 

38 


78 

66 

222 

185 

208 

151 


-201 

-187 

-120 

43 

-98 

-186 


335 

307 

100 

-64 

112 

191 


10i4 

1058 

852 

600 

690 

729 

(Pz.I’e) 

555 

571 

512 

425 

460 

363 


79 

91 

11 

-8 

19 

65 
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Thirty-day rr.ean spectral distributions of K(n) / P(n), 
M'li) , and C(n) for the control and anomaly runs for tl 
first month in the zonal section 2 8 U-GO'^N. (A) Spec- 
trv'il component of eddy kinet.ic energy, K(n) . (B) Spec- 
tral componenc of ecidv avjiilable potential energy, P(n) 
(C) conversion from zonal kinetic energy to 7.{vi) , M(a) 
(B) Co:ivcr.vioi\ from P(:'i) to K(n', C(n). 
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and C(n) for each month of the control and anomaly runs in 
the same zonal section. K(n) and P(n) are the spectral com- 
ponents of K and ? , respectively , ' and K(n) and C(n) the 

spectral components of ^ ) ' respectively. 

Table 3 Indicates that In the first month the 
effects -in the northern middle latitudes are fairly small 
compared with the other two months. In the first month, the 
trainslent SST anomaly appears to cause an increase in Pg and 
Gg , as was expected, and a slight reduction in i together 
with slightly increased conversions from to Pg , Pg to Kg , 
and Kg to K^. However, during the first month the spectral 
energetics in the northern middle latitudes are greatly dis- 
turbed by the SST anomaly (Pig. 17). As might have been 
expected, the anomaly appears to cause an increase in the 
eddy activity of the intermediate scales. In the intermediate 
scales (n=4“10), eddy energies, K(n) and P(n), and eddy con- 
version, C(n), are increased in the anomaly run (Figs. 1?A, 
17B, 17D), On the northern side of the warm pool, the sea 
surface temperature gradients are increased in the first 
month by the prescribed SST anomaly. This increases the 
barocllnlc instability and eddy conversion around the inter- 
mediate scales. The more active cyclone waves then transfer 
more kinetic energy to zonal flow to maintain the zonal wind 
system (Fig. 1?C). Although not shown in the figures, the 
spectral analysis reveals that the increase of Gg takes place 
mainly in wave number 1 and 2. This increase is due primarily 
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to the increased latent heat released by condensation. 

Another noticeable feature occurlng in the first 
month is that the eddy activity of the ultra-long waves is 
markedly reduced in the anomaly run, especially wave number 
2 (Fig, 17). Because most of the increase of eddy activity 
around the intermediate scales is compensated by the reduc- 
tion In ultra-long waves, the net eddy activity shows only 
a minor Increase In the northern middle latitudes (Table 3). 
Heat sources and slTiks of planetary dimensions are very 
important for determining large-scale flow patterns (e,g,, 
Smagorinsky, 1953). Large-scale flow patterns are generally 
dominated by the long planetary waves. The SST anomaly in 
the first month probably cause a considerable change in the 
distribution of large-scale heat sources and sinks. This 
then produces a larger dynamical response in the ultra-long 
waves. However, the way that ultra-long waves respond to 
SST anomalies is not so easy to anticipate, unlike the res- 
ponses of the Intermediate cyclone waves, which can be in- 
ferred from barocllnlc Instability theory. 

In the second month, as indicated In Table the 

responses in the northern middle latitudes are more apparent. 

The anomaly causes eddy generation, G , to reduce to zero 

£ 

In the anomaly run. It also causes ' 

and Wg to change sign. The change of < indicates 

a direct energy conversion replaced the normal indirect energy 
conversion. The reduction In the strength of the upward 



67 


branch of the Perrel cell Is responsible for this change 

(Fig. 8A). The reversal in sign of .( K„,K ^ indicates an 

' r. z ' 

anomalous conversion from zonal to eddy kinetic energy. 

These changes together with a reduction in Kg , Pg , { , 

and ( P_,P„ ) clearly indicates that the transient SST anomaly 
^ E ' 

has a effect of reducing eddy activity in the northern middle 
latitudes after a period of one month. The spectral analysis 
reveals that the reduced eddy activity takes place at all 
wave numbei's, with the major reduction in the longer waves 
(Pig. l8). Prom Fig. 18 it can be seen that K(n), P(n), and : 
C(n) (also eddy generation, G(n), but not shown In the 
figures) decrease over almost all waves, with the dominant 
reduction in the longer waves. In a barotroplc atmosphere, 
the unstable eddy receives its energy from zonal kinetic 
energy (e.g,, Kuo, 19^9; Brown, 1969). Fig. 18C illustrates 
the fact that the anomaly causes the eddy of wave number 2 
to become less barotropically unstable, but causes the eddies 
of wave number l and 3 to become more barotropically unstable. 
The reversal in sign of ( primarily to the 

changes in wave number l and 3, 

In the third month, as can be seen from Table 3, 

^ northern middle latitudes are 

reduced in the anomaly run, but , Kg , and ^ K^,K^ ) are 
increased by the transient Northern Pacific warm pool. The 
anomaly also intensifies the indirect Perrel cell and causes 
eddy pressure interaction to transport more energy cut of 



68 


northern middle latitudes. During the third month, the eddy 
activity is generally increased in the northern middle lati- 
tudes. In the wave-number domain, this increase takes place 
from n^2 to n=7 (Fig. 19). Fig. 19 illustrates an increase 
of K(n), P(n), and C(n), and a decrease of H(n) at these 
waves. .Among those waves, eddy conversions at wave number 2 
and 6 are particularly Increased In the anomaly run (Fig. 19D). 
Also, all the waves were found to transfer kinetic energy to 
zonal flow in the anomaly run, while zonal flow supplies 
energy to eddies of wave number 1, 3» 5 £oid 6 in the control 
run (Fig. 190). This is why <Kg,K^) Is increased in the 
anomaly run. 

(B) Effects in the southern middle latitudes 


The thirty-day mean energetics for each month of 
the control and anomaly runs in the zonal section 28°S - 60:°S 
are shewn in Table 4. The thirty-day mean spectral distribu- 
tions of K(n), P(n), K(n), and C(n) for each month of the 
control and anomaly runs in the same region are illustrated 
in Figs. 20, 2i, and 22, The gross regional energetics both 
in the spatial and wave-number domains shovf that the effects 
in the southern middle latitudes are very small during the 
first month (Table Fig, 20). 

In the second month, the most noteworthy effects 
of the transient SST anomaly in the southern middle latitudes 
are the reversal in alpsn of , and the large decrease 

in wV (Table 4). The anomaly also inc7‘-easer, the conversioris 
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Table Thirty-day mean energetics for each month of the 

control and anomaly runs In the zonal section 

28^S - 60^S. The unit of energy content is 10^ 
—2 

Joules m and that of energy transformation is 

•2 —2 
10 watt m • 


Month No. 

Control 

1 

Anomaly 

1 

Control 

2 

Anomaly 

2 

Control 

3 

Anomaly 

3 


209 

226 

286 

305 

250 

259 


84 

80 

89 

98 

78 

47 

^z 

153 

152 

l68 

180 

l6i 

168 


73 

72 

88 

94 

85 

57 

®z 

119 

127 

210 

199 

185 

174 


53 

46 

31 

13 

35 

28 

«z 

l4l 

l40 

202 

219 

183 

ll^5 


-27 

-85 

l60 

34 

6i 

-ll6 

( ) 

101 

91 

95 

138 

i49 

163 

(Pe.Ke) 

665 

644 

556 

573 

628 

56l 


99 

108 

78 

-38 

129 

87 

(Vz ) 

47 

46 

36 

59 

51 

82 
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Same as Fig. 20 , for the third month. 
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{ ^z'^z / ( ^E*^z / ' eddy enercries in the anomaly 

run. Fig. 21 illustrates that ultra-long and intermediate 
■waves are fairly disturbed by the SST anomaly. The anomaly 
appears to increase particularly the eddy activity of wave 
number 5- For* example, K(5), P(5)> and C(5) are markedly in- 
creased, and H(5) markedly decreased in the anomaly run. 
Although not shown in the figures, the eddy of wave number 5 
was found to transfer more available potential energy to 
zonal in the anomaly run. This Increase contributes much 
to the reversal in sign of ( Fz»^E ) southern middle.;, 

latitudes. 

In the third month, the most noteworthy responses 
td the anomaly in the southern middle latitudes are the re- 
versal in sign of Wg and the large decrease in the eddy 
energies (Table 4), The anomaly causes Wg in the anomaly 
run to transport energy out of the southern middle latitudes 
Instead of importing energy into the southern middle lati- 
tudes as in the control run. The conversions / P_,K \ and 

\ E / 

{ Pg.Pg ) are also decreased by the transient warm pool. 

Figs, 22k and 22B illustrate that eddy energies decrease 
over all waves in the anomaly run, . The reduction of eddy 
conversion also occurs over most of the waves (Fig. 22D). 

In the anomaly run all waves in the southern middle lati- 
tudes were found to be barotroplcally stable, while the eddy 
of wave number 3 the control run is not (Fig. 22C). 

In general, the analysis of the regional energecics 



both In the spatial and wave-number domains show that the 
middle latitudes of both hemispheres can indeed exhibit a 
great response to a transient North Pacific SST anomaly. 
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8, Conclusions and further remarks 

Two aspects of the problem under investigation are 
(l) the ability of the two-level Mintz-Arakawa model to simu- 
late the atmospheric energetics, and (2) the possible effects 
of a transient North Pacific SST anomaly on the energetics of 
the atmosphere. For the first problem, the energetics of the 
control run, which was designed to simulate normal climate, 
were compared with the observations. For the second problem, 
the energetics of the anomaly run, corresponding to a tran- 
sient anomalous SST patteim, were compared with the control 
run. 

Except in the tropics, the two-level Mintz-Arakawa 

model exhibits a reasonable and realistic energy budget, 

o 

For the area north of 20 N, except for the dlabatic genera- 
tion of eddy available potential energy, the energy budgets 
both in the spatial and wave-number domains show good agree- 
ment between the control rim and the observations. In the 
tropics, the convective precipitation was found to have a 
dominating Influence on the model energetics. 

One major defect of the model simulation is caused 
by the parameterization of convection, which produces too 
much convective precipitation in the tropics. The excessive 
rainfall is accompanied by the release of too much latent 
heat of condensation, and hence magnifies the diabatic genera- 
tion of available potential energy in the tropics. It also 
increases the vertical velocities v;hich in turn intensify 
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the conversion from available potential to kinetic energry. 

These results further affect other atmospheric energy processes 
in the tropics. Thus, improvement on the parameterization 
of convection may bring the simulated tropical energetics 
closer to the reality. 

The analysis of the spectral energy budget of the 
control run for the area north of 20*^N reveals that the 
spectral components of eddy energies, P(n) and K(n), are 
generally smaller in the model than in nature, as are the 
conversions from P_ to P(n). The coarse horizontal resolu- 
tion in the model might be responsible for these shortcomings. 
The model was found to produce excessively high temperature 
In the tropics, prssumely due to the excessive convective 
precipitation in the tropics. This bring about a much larger 
meridional temperature gradient in middle latitudes in the 
model. Therefore, the fact that the conversions from to 
P(n) are too small in the model implies that the eddies are 
less effective in transporting sensible heat poleward. 

In the SST anomaly experiment, the imposed positive 
anomalous pattern has tv;o immediate effects in the first month. 
One is to augment the sensible heat transfer from sea to air 
and t'na evaporation over the local anomaly region. This in 
turn increases the released latent heat and hence increases 
the diabatic generation of eddy available potential energy, 

Og. The other is to increase the sea surface temperature 
gradients, especially on the northern side of the warm pool. 
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This then- increases the baroclinlc instability and eddy con- 
version, C(n), around the intermediate scales. Since the 
mean motion of the atmosphere is mainly parallel to the lati- 
tude circles, and the cyclone paths are toward northeast, 
these effects spread to the other parts of the northern 
middle latitudes. These two effects further intensify other 
eddy activities, e.g,, K(n), P(n), Q(n), and K(n), around the 
intermediate scales in the northern middle latitudes (Fig. 17). 

However, an unexpected result in the first month 
is that the anomaly appears to cause a marked reduction in 
the eddy activity of the ultra-long waves, especially wave 
number 2, in the northern middle latitudes (Fig. 17). Be- 
cause most of the increase of eddy activity around the 
intermediate scales Is compensated by the reduction in 
ultra-long waves, the net eddy activity shows only a minor 
Increase in the northern middle latitudes (Table 3)* 

Although the first month's gross regional energetics in the 
spatial domain are slightly affected by the anomaly (see 
Tables 3 ^)i the spectral energetics in the northern 

middle latitudes appear to be considerably disturbed by the 
SST 6inoraaly. 

In the second month, with the removal of the 

anomalous SST pattern, the poi-tion of the atmosphere above 

the first month's warm pool tends to experience a reduction 

of surface sensible heat exchuunge and evaporation. This 
reduces the condensation and hence decreases Gp. 


It also 
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decreases the vertical velocities, end this in turn reduces 
the eddy conversion ^ These effects further affect 

other energetics and induce significantly large response 
in the northern middle latitudes. As shown In Table 3» eddy 
energies, Pg and K^, , and the conversion { also 

reduced in the anomaly run, and the conversions and 

^ ) , and the eddy pressure interaction, Wg , in the 

anomaly run are changed In sign as compared with the control 
run. The reduction in weaker Indirect cell 

further imply that the poleward eddy transports of sensible 
heat and momentum may be reduced by the anomaly. In general, 
the transient SST anomaly causes a reduction in eddy activity 
in the northern middle latitudes. The reduced eddy activity 
was found to take place at all wave numbers, with the ma^or 
reduction In the longer waves (Pig. l8). 

The southern middle latitudes in the second month 
are also considerably disturbed by the anomaly, e,g,, a 
reversal in sign of < Pz^Pg) » and an increase in < * 

( ) , Pg , and Kg (Table 4), Longer waves are fairly 

disturbed by transient warm pool, especially wave number 5 
showing a marked increase in eddy activity (Pig. 2l). 

The global energetics in the third month are still 
fairly disturbed by the transient SST anomaly, Acco 3 ?ding to 
Tables 3 and 4 and Figs. 19 and 22, the eddy activity gen- 
erally increases in northern middle latitudes tut decreases 
in southern middle latitudes, and the variations of eddy 
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activities take place primarily in the longer waves. 

The third month's tropical energetics were found 
to be severely distorted in the aiionialy run compared with 
the control run (Figs. 5, 6, 12 - l6). This is due primarily 
to the extremely large reduction of precipitation at 2°S. 

Due to the unsatisfactory parameterization of convection in 
the model, this tropical effect Is not considered as a true 
reflection of nature, but rather as a computational artifact 
caused by a model defect. In order to simulate more reliably 
the effects of SST anomalies in the tropics, an improvement 
in the treatment of convection is badly needed. (Note that 
this improvement was included in Arakawa (1972) 3-l&y®r model, 
and has been retained in GISS 9-layer model (Somerville et al., 
1'974)), 

Although the anomalous SST pattern persists for 
only one month, the results indicate that the transient SST 
variations of reasonable magnitude in the North Pacific Ocean 
can induce a disturbing effect on the global energetics both 
in the spatial and wave-number domains. This implies that 
modest variations in the sea surface temperature may cause a 
significant change in the monthly or seasonal synoptic patterns 
over the whole earth, and that a precise knowledge of the 
time-dependent ses surface temperature field may be important 
for the long-range weather prediction. It remains to be seen 
v^'hethsr this is a property only of the model, or if it is 
indeed a true reflection of nature. 
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